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Abstract. We apply the subtracted kernel method (SKM), a renormalization 
■ approach based on recursive multiple subtractions performed in the kernel of the 

scattering equation, to the chiral nucleon-nucleon (NN) interactions up to next-to- 
next-to-leading-order (NNLO). We evaluate the phase-shifts in the 1 So channel 
at each order in Weinberg's power counting scheme and in a modified power 
counting scheme which yields a systematic power-law improvement. We also explicitly 
" demonstrate that the SKM procedure is renormalization group invariant under the 

£*"»« . change of the subtraction scale through a non-relativistic Callan-Symanzik flow 

ov equation for the evolution of the renormalized NN interactions. 

1. Introduction. 

The issue of the non-perturbative renormalization of nucleon-nucleon (NN) interactions 
in chiral effective field theory (ChEFT) has been intensively investigated by many 
authors (for detailed reviews, see e.g. Refs. [H 121 El H]), generating a great deal of 
discussion and debate regarding the consistency of the approach originally proposed by 
Weinberg [3 El [7]. The standard procedure for the non-perturbative renormalization 
of NN interactions in the context of Weinberg's approach to ChEFT can be divided in 
two steps [2]. In the first step, one has to solve a regularized Lippmann-Schwinger (LS) 
equation for the scattering amplitude by iterating the effective NN potential truncated 
at a given order in the chiral expansion, which includes long-range contributions from 
pion exchange interactions and short-range contributions parametrized by nucleon 
contact interactions. The most common scheme used to regularize the ultraviolet 
(UV) divergences in the LS equation is to introduce a sharp or smooth momentum 
cutoff regularizing function [21 [3] that suppresses the contributions from the potential 
matrix elements for momenta larger than a given momentum cutoff scale (multi- 
pion exchange interactions also involve UV divergent loop integrals which must be 
consistently regularized and renormalized). In the second step, one has to determine 
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the strengths of the contact interactions, the so called low-energy constants (LEC's), 
by fitting a set of low-energy scattering data. Once the LEC's are fixed at a given 
momentum cutoff scale, the LS equation can be solved to evaluate other observables. 
The NN interactions can be considered properly renormalized when the predicted 
observables are (approximately) independent of the momentum cutoff scale within the 
range of validity of ChEFT [31 [H] - 

The state of the art chiral NN potentials available to date, constructed within the 
framework of Weinberg's approach to ChEFT, are the next-to-next-to- next-to-leading- 
order (N 3 LO) potentials of Epelbaum, Glockle and Meissner |9] and of Entem and 
Machleidt [10J. Both these potentials provide a very accurate description of NN 
scattering data below laboratory energies -Elab ~ 300 MeV, with a x 2 /d.o.f ~ f 
comparable to that obtained by high-precision phenomenological potentials such as 
the Nijmegen II [H] and the Argonne V18 JL2], and have been successfully applied 
in many nuclear structure and reaction calculations. Furthermore, the leading chiral 
three-nucleon (3N) forces at NNLO have been derived [T3l [H] and applied with fair 
success in calculations of few-nucleon reactions, light and medium nuclei, and infinite 
nuclear and neutron matter. Subleading chiral 3N forces at N 3 LO [T5l [TBI [TT] and 
N A LO [HI [19j [20] have been recently worked out and are expected to resolve some 
challenging nuclear structure and reaction problems that remain open. For a review 
and a comprehensive list of references on the applications of ChEFT to few- and many- 
nucleon systems, see e.g. [21 El II]- 

In spite of its remarkable phenomenological success, the theoretical basis of 
Weinberg's approach has been criticized. In particular, conceptual questions have been 
raised regarding the formal inconsistency between Weinberg's power counting scheme 
(based on naive dimensional analysis) and the non-perturbative renormalization of the 
NN interactions [21 E]- Such inconsistency arises from the non-renormalizability of 
ChEFT, which is a consequence of the highly singular nature of the interactions in 
the chiral expansion of the NN potential at short distances. The iteration of the 
NN potential truncated at a given order in the chiral expansion (using the LS or 
the Schrodinger equation) generates higher-order UV divergences, and hence strong 
regularization scale dependencies, which cannot be absorbed by the contact interactions 
(counterterms) included in the potential at that same order. Thus, upon momentum 
cutoff regularization, the limit of infinite cutoff cannot be taken while keeping the 
amplitudes finite and cutoff- independent, i.e. renormalization group invariant. So far, 
the issue of proper non-perturbative renormalization of chiral NN interactions remains 
controversial. From one side, in the successful ChEFT calculations based on Weinberg's 
approach, implemented by using the chiral N 3 LO potentials mentioned above, finite 
cutoffs are typically chosen in the range ~ 450 — 600 MeV and nearly cutoff- independent 
stable results are obtained provided the cutoff is varied only over a narrow window 
[21 El II]- A justification for such a setup with basis on Lepage's view of renormalization 
in cutoff EFT [8] has been presented in the works of Epelbaum and Meissner [21] and 
Epelbaum and Gegelia [22]- Still within the framework of finite cutoff calculations, it 
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has been shown that the explicit inclusion of A — isobar degrees of freedom in ChEFT 
improves the convergence of the chiral expansion for the NN interactions as compared 
to the theory with only pion and nucleon degrees of freedom [231 EU EE, EHl ETJ EH] . From 
the other side, the criticism to the consistency of Weinberg's approach, and in particular 
to the narrow range of cutoffs for which it can provide renormalization group invariant 
results, led to the proposal of alternative renormalization approaches and power counting 
schemes such as those by Kaplan, Savage an Wise [291 EHl EI] , van Kolck [32] , Nogga 
et al. [33], Birse PH EHl EE] , Valderrama (HTJEHIES], Beane et al. gQJ SI] and Long 
et al. [121 H3j HI]. In this scenario, methods based on Wilson's renormalization group 
[4"5l SE1 H7] have been successfully applied to analyze the scale dependence of chiral 
NN potentials, both in momentum space [3H [351 EE1 SHI SI] and in coordinate space 
[50] I5T1 l52l l53l 154] , providing a better understanding of the interplay between the power 
counting and the non-perturbative renormalization in ChEFT. Although much progress 
has been made in this direction, the construction of a consistent framework for the non- 
perturbative renormalization of NN interactions still remains an open problem [55] . 

An alternative approach to the non-perturbative renormalization of NN 
interactions that has also been explored is the subtracted kernel method (SKM) 
[56| [571 1581 1591 1601 loll 1621 1631 164"! [65] in which, instead of using a momentum cutoff 
regularizing function, the LS equation is regularized by performing multiple subtractions 
in the kernel at a given energy scale. A closely related subtractive renormalization 
approach is described in Refs. [66[ [67[ [68] , although in that case a sharp momentum 
cutoff is also introduced. In this work we apply the SKM approach to the chiral NN 
interactions up to next-to-next-to- leading-order (NNLO). We consider the scattering 
of two nucleons in the 1 5'o channel and analyze the errors in the phase-shifts calculated 
at each order in Weinberg's power counting scheme (WPC) and in a modified power 
counting scheme (MPC) based on the promotion of contact interactions which yields a 
systematic order- by-order power-law improvement. We also show, by explicit numerical 
calculations, that the SKM procedure is renormalization group invariant under the 
change of the subtraction scale, provided the renormalized interactions are evolved 
through a non-relativistic Callan-Symanzik (NRCS) flow equation [69l [70] . 

2. SKM approach for the NN system. 

We start by considering the chiral expansion for the effective NN potential in Weinberg's 
power counting scheme (WPC) [21 [3]. In a partial-wave relative momentum space basis, 
the matrix-elements of the NN potential in the 1 S'o channel up to N 3 LO are given by 




(p,pO 

°(P,P') 

TL °(P,P') 

LO (p,p') 



+ V 3 f LO (p,p') + CJ p 2 p' 2 + C^+p' 4 ) , (4) 
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where the coefficients Cj stand for the strengths of the contact interactions and an 
obvious notation is used for the pion-exchange interactions. 

The matrix-elements of the leading-order (LO) OPE potential in the 1 Sq channel 
are given by 

32iij* \ J-i p l + p' 2 — 2pp'x + m 2 J 

where g a , f n and m n denote, respectively, the axial coupling constant, the pion weak- 
decay constant and the pion mass. The higher-order OPE terms include corrections from 
pion loops and counter term insertions, which only contribute to the renormalization 
of coupling constants and masses. In this work, we use g a = 1.25, f n = 93 MeV and 
= 138 MeV. The two-pion-exchange (TPE) potential is taken from Ref. [9]. 
Consider the formal LS equation for the T-matrix of a two-nucleon system, which 
can be written in operator form as 

T(E) = V + V Gq(E) T(E) , (6) 

where E is the energy of the two-nucleon system in the center-of-mass frame, V is the 
effective NN potential and Gq(E) = (E — Hq + ie)~ x is the free Green's function for 
the two-nucleon system with outgoing-wave boundary conditions, given in terms of the 
free hamiltonian Hq. Both pion-exchange and contact interaction terms can lead to UV 
divergences when the effective AW potential V at a given order in the chiral expansion 
is iterated in the LS equation, requiring a regularization and renormalization procedure 
in order to obtain well-defined finite solutions. 

In the standard cutoff renormalization scheme the formal LS equation, Eq. (JSJ), 
is regularized by multiplying the effective iViV potential V with a momentum cutoff 
regularizing function. The common choice is an exponential f(p) = exp[— (p/A) 2r )] 
(with r = 1, 2, . . .), where A is a cutoff parameter, such that 

V(p,p') -> V A (p,p') = exp[-(p/A) 2r )] V(p,p') exp[-(p'/A) 2r )] . (7) 

In the SKM approach, a regularized and renormalized LS equation for the T- 
matrix is computed through an iterative procedure which involves recursive multiple 
subtractions performed in kernel at a given energy scale. For a general number of 
subtractions n, we define a n-fold subtracted LS equation given in operator form by 

T (n) {E) = yin) {E) + y(n) {E) _ ^ ^n) {E) ^ (g) 

where fx is the subtraction scale, V^ n \E) is called "driving term" and G+(E; — /i 2 ) is 
the n-fold subtracted Green's function, defined by 

G+(E; -u 2 ) = (^4-) G+{E) = F n (E; -u 2 ) G+(E) . (9) 



Note that we choose a negative energy subtraction point — « 2 , such that the free Green's 
function Gq (— u 2 ) is real. 

The n-fold subtracted LS equation, Eq. (jSJ), has the same operator structure as 
the formal LS equation, Eq. (jBJ), but with the effective AW potential V replaced by 
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the driving term V^^E) and the free Green's function Gq(E) replaced by the n-fold 
subtracted Green's function G^(E\ —fi 2 )- Moreover, regularization is achieved not by 
using a momentum cutoff regularizing function, but instead through the form factor 
E n (E; —fi 2 ) introduced by the n-fold subtracted Green's function, which is built up 
as part of the iterative procedure. The driving term encodes the physical information 
apparently lost due to the removal of the propagation through intermediate states at the 
subtraction scale \i. Thus, once the driving term is determined (by fixing the strengths of 
the contact interactions at the subtraction scale \i) the subtracted LS equation provides 
a renormalized finite solution for the T-matrix at any given energy E. 

The driving term Vjf 1 ' (E) is recursively constructed through an iterative procedure 
starting from V^ 1 \E = — fi 2 ) = T^(E = —/J 2 ), which is replaced by an ansatz for the 
T-matrix at the subtraction scale, T(—fi 2 ). The recursion formula (with m = l,...,n) 
is given by 

V^(E) = V^(E)+V^ ins , (10) 

where 

Vl 1 m \E) = [l - (-fi 2 - E) m ~W^- 1 \E)Gt(-fi 2 ) m ]'W^- 1) (E) , (11) 

and the term V^ ing (E) contains the higher-order singular interactions that generate 
divergent integrals which can be regularized by performing m subtractions. One should 
note from Eq. fflUj) the that the driving term V^ m \E) at each iteration is derived in 
two steps. First, we calculate V^ m \E) from V^ m-1 '(-E), solving an integral equation 
obtained by manipulating Eq. ({TTI) : 

Vjjri{E) = V^ m - 1) {E) + {-fi 2 -E) m - 1 V^ m - l \E) G+{-fi 2 ) rn V^ m \E) .(12) 

Then, we introduce the corresponding higher-order singular interactions in the driving 
term by adding V^ ing (E). 

For convenience, in this work we implement the SKM procedure using the i^-matrix 
instead of the T-matrix. The LS equation for the i^-matrix in the 1 5'o channel with n 
subtractions is given by (here and in what follows we use units such that h = c = M = 1, 
where M is the nucleon mass) 



K^(p,p'-k 2 ) = V^\ P , P '-k 



(is) 

where k = y/E is the on-shell momentum in the center-of-mass frame and V denotes the 
principal value. Note that the n-fold subtracted Green's function introduces an energy- 
and /i-dependent form factor proportional to q~ 2n in the kernel of the subtracted LS 
equation which regularizes UV power divergences up to order (9(g 2n_1 ), thus effectively 
acting like a smooth momentum cutoff regularizing function. 

At each order in the chiral expansion we perform the minimum number of 
subtractions n necessary to render a finite solution for the subtracted i^-matrix. The 
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corresponding driving term is computed through the iterative procedure described 
above, starting from the ansatz for the LO driving term V^ 1 ^ given by 

V™M; -/x 2 ) = V™(p,p') + CM . (14) 

One should note that only one subtraction is enough to get a finite result for the K- 
matrix at LO, since when iterated in the LS equation the LO potential generates UV 
power divergences of order 0(q l ). 

At next-to-leading-order (NLO) we have to construct a 3-fold subtracted kernel LS 
equation in order to get a finite result for the L'-matrix, since when iterated in the LS 
equation the NLO potential generates UV power divergences up to order 0(q 5 ), even 
thought it contains only terms up to order 0(q 2 ). To obtain the NLO driving term, 
Vj?\ we first calculate from V^ l \ 

vf> (p, p'; k 2 ) = vm ( P , p'; k 2 ) - i jf d q e vm (p , q - k 2 ) 

Then, we calculate from V^ 2 \ 

V«>{p,p'-,k 2 ) = Vf\p,p\k 2 ) -\[dq q 2 V^(p,q;k 2 ) 

x J^W 9 ^ {q ' p '^ h (16) 

and add the NLO interactions: 

Vf\p, P '; k 2 ) = V®(p,p'; k 2 ) + V 2 N n LO (p,p') + CM {p 2 + p' 2 ) . (17) 

At NNLO we have to construct a 4-fold subtracted kernel LS equation in order to 
get a finite result for the i^-matrix. Even thought the NNLO potential contains only 
terms up to order 0(q 3 ), when iterated in the LS equation it generates UV power 
divergences up to order 0(q 7 ). To obtain the NNLO driving term, V^ 4 \ we first 
calculate Vffi from 



9 roo 

V^(p,p';k 2 ) = if >(p,p';k 2 ) -± jf dq q 2 V«>(p,q; k 2 < 

x j^3^ 4, («.p';* 2 ). (18) 

and then we add the NNLO interactions: 

V^(p jP >; k 2 ) = ?W(p,j/; k 2 ) + V» NLO (p,p') . (19) 

The renormalized strengths Cj(/i) of the contact interactions included in the driving 
term at each order in the chiral expansion are fixed at the subtraction scale jx 
by fitting data for low-energy scattering observables, thus encoding the input physical 
information. Instead of the usual matching of scattering data at discrete values of 
the on-shell momentum k, we follow the procedure described by Steele and Furnstahl 
[7T1 172] . which is numerically much more robust. Here, we use as "data" the values of 
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the inverse on-shell i^-matrix evaluated from the solution of the LS equation with the 
Nijmegen-II potential [11] for a spread of very small momenta k (< 0.1 fm _1 ) and fit 
the difference between such data and the inverse on-shell i^-matrix evaluated from the 
solution of the n-fold subtracted LS equation (Eq. ({TBI ) with the driving term to 
an interpolating polynomial in fc 2 //i 2 to highest possible degree, 

A(l/K) = 1/K mj (k, k; k 2 ) - 1/K^(k, k; k 2 ) 
k 2 k 4 

= A + A 2 — + A A — + ... . (20) 

The coefficients Aj are then minimized with respect to the variations in the renormalized 
strengths Cj(/i). One should note that the procedure outlined above is equivalent to fix 
the renormalized strengths Ci(fi) by fitting the experimental values of the parameters 
of the effective range expansion (ERE) [73] to a given order in k 2 , since at very low 
energies the phase-shifts provided by the ERE agree very well with those obtained from 
the Nijmegen-II potential, i.e. 

- l/K NU (k,k;k 2 ) = k cot6 mj (k) ~ k cot 5 ERE (k) 

= - ~ + \r e k 2 + v 2 k* + v 3 k 6 + . . . , (21) 
where a is the scattering length, r e is the effective range and Vi are the shape parameters. 

3. Power counting for the NN system in the 1 So channel. 

Once the renormalized strengths Cj(/x) are fixed at the subtraction scale /z, and so the 
driving term is determined, we can calculate the NN scattering observables for any 
given energy from the numerical solution of the n-fold subtracted LS equation for the 
fT-matrix (Eq. ( fl3l) ). In log- log plots for the relative errors in the observables ("Lepage 
plots"), we expect to obtain straight lines with slopes given by the dominant power 
of k 2 //j, 2 in the errors. A systematic power-law improvement in the predictions of the 
observables is expected as more contact interactions are included in the driving term 
V( n > and the corresponding renormalized strengths Cj(u) are fixed. The impact of the 
included long-range pion-exchange interactions on the scaling of the errors is, however, 
much less transparent and so must be verified through the explicit numerical calculation 
of the observables. 

In the left panel of Fig. [1] we show the phase-shifts in the 1 So channel as a function 
of the laboratory energy (-Elab) calculated from the numerical solution of the n-fold 
subtracted LS equation for the K- matrix at LO (n = 1), NLO (n = 3) and NNLO 
(n = 4), with the renormalized strengths Cj(/i) of the contact interactions included at 
each order fixed at the subtraction scale fi = 1.03 fm _1 . In the right panel we show the 
log-log plots for the corresponding relative errors with respect to the Nijmegen-II phase- 
shifts. For the phase-shifts calculated at LO, where the driving term consists of the OPE 
interaction plus the non-derivative contact interaction, we obtain the well-known strong 
deviation from the Nijmegen-II results, with the relative errors scaling like 0(k 2 //j, 2 ) at 
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Figure 1. (Color on-line) Phase-shifts in the 1 S'o channel calculated from the solution 
of the subtracted LS equation for the if-matrix at LO, NLO and NNLO both in 
Weinberg's power counting scheme (WPC) and in the modified power counting scheme 
(MPC) with the renormalized strengths Ci(/i) fixed at fj, = 1.03 fm _1 (left) and the 
corresponding relative errors with respect to the Nijmegen-II phase-shifts (right). 



low-energies. For the phase-shifts calculated at NLO, where a TPE interaction term 
and the second-order derivative contact interaction are included, we obtain the expected 
power-law improvement, with the relative errors scaling like 0(k A /n 4 ) at low-energies. 
At NNLO in the WPC scheme, where only a TPE interaction term is included but 
no new contact interaction, we obtain a better overall agreement with the Nijmegen-II 
phase-shifts but there is no power- law improvement. 

We then consider a modified power counting scheme (MPC) in which the fourth- 
order 1 S'o channel derivative contact interaction C' 4 (fi) p 2 p' 2 is promoted from N 3 LO 
to NNLO and, therefore, added to the driving term Vjf> given by Eq. (|T9|) . As one 
can see from the error plots shown in the right panel of Fig. [U at NNLO in the MPC 
scheme the relative errors in the phase-shifts scale like 0(k 6 /fi 6 ) at low-energies and so 
we obtain a systematic order- by-order power- law improvement. 

Modifications of the WPC scheme based on the promotion of contact interactions 
have been considered in several works jSSlEllESlESlEZlESlESlHSlllSISl- In the 
power counting scheme proposed by Nogga, Timmermans and van Kolck (NTvK) |33j, 
which is also used in the nuclear matter calculations described in Ref. [71], higher- 
order contact interactions are promoted to LO in low angular momentum partial-wave 
channels where the OPE tensor interaction is singular and attractive (namely, in the 
3 Pq, 3 P2 — 3 F 2 and 3 D 2 channels), such that reasonably cutoff independent results are 
obtained in LO calculations for cutoffs varying in the range 2 ... 20 fm _1 . A criticism 
of the NTvK scheme was made in Ref. [21] where, based on an analysis of several NN 
scattering observables calculated at LO both in the NTvK and the WPC scheme, it has 
been argued that the use of larger cutoffs and the modifications of the WPC scheme as 
proposed in Ref. [33] do not improve the results and can even lead to discrepancies in 
certain partial-wave channels. 
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Our MPC scheme resembles that implemented by Valderrama [371 ESI ES] within 
the framework of a perturbative treatment of the chiral TPE interactions at NLO 
and NNLO, which can be regarded as an extension of the NTvK scheme to sub leading 
orders. In Valderrama's approach, the LO phase-shifts are calculated non-perturbatively 
from the solution of the Schrodinger equation in coordinate space with the LO potential 
(OPE interaction plus the non-derivative contact interactions) iterated to all orders. 
The contributions to the phase-shifts from the subleading order terms of the potential 
are included perturbatively within the distorted- wave Born approximation (DWBA), 
resulting in a set of power counting rules (based on the requirement of renormalizability 
of the scattering amplitude) which are in agreement with the modifications of the WPC 
scheme determined in the renormalization group analysis made by Birse [34"] , 135 } 136] (with 
some minor differences). A perturbative treatment of the subleading order interactions 
is also considered in Refs. 0QHHlg2lH3lH3lCSlCSllZZllZB]. Here, we consider 
the application of the SKM approach to the chiral NN interactions within a non- 
perturbative renormalization framework and so we calculate the phase-shifts by iterating 
the full NN potential (LO plus subleading order terms) to all orders in the subtracted 
LS equation. One should also note that in Valderrama's approach the power counting 
scheme is such that it requires the promotion of the 1 Sq channel contact interaction 
C4 (p 4 + p' 4 ) from N 3 LO to NLO, while in our MPC scheme we promote the 1 S'o 
channel contact interaction C' 4 p 2 p' 2 from N 3 LO to NNLO. 

In principle, we could promote either of the two 1 5'o channel fourth-order derivative 
contact interactions available at N 3 LO in the WPC scheme. We choose to promote the 
C' A (fi) p 2 p' 2 contact interaction because when iterated in the standard LS equation it 
generates additional UV power divergences only up to order 0(q 6 ) (due to the presence 
of the NNLO TPE interaction term), which can be regularized by performing the 
same four subtractions as required for the NNLO calculation in the WPC scheme. 
The iteration of the C^fi) (p 4 + p' 4 ) contact interaction, on the other hand, generates 
additional UV power divergences up to order 0(q 9 ), and so would require at least five 
subtractions. Moreover, with the C^/i) {p 4 + p' A ) contact interaction left to be included 
at N 3 LO, there will be no need to promote a higher-order contact interaction to obtain 
the power-law improvement at this order. In this way, using such a prescription we 
would obtain in our MPC scheme the same number of 4 Sq channel contact interactions 
at N 3 LO as predicted in the WPC scheme, which is the one employed in the construction 
of the state of the art chiral N 3 LO chiral potentials described in Refs. |9j and [10J. 

In the left panel of Fig. [2] we show the phase-shifts in the 1 S'o channel as a function 
of £"lab calculated from the solution of the subtracted LS equation for the .fT-matrix 
at NNLO in the MPC scheme, with the renormalized strengths Cj(/i) (i = 0, 2, 4) 
fixed at several values of the subtraction scale /i. In the right panel we show the log-log 
plots for the corresponding relative errors with respect to the Nijmegen-II phase-shifts. 
As one can see, the range of energies for which the SKM procedure provides a good 
description of the phase-shifts increases with fi. For fi = 1.03 fm _1 , the SKM results 
are in good agreement with the Nijmegen-II phase-shifts up to £"lab ~ 60 MeV. 
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Figure 2. (Color on-line) Phase-shifts in the 1 Sq channel calculated from the solution 
of the subtracted LS equation for the if-matrix at NNLO in the modified power 
counting scheme (MPC) for several values of the subtraction scale fj, (left) and the 
corresponding relative errors with respect to the Nijmegen-II phase-shifts (right). 

It is important to observe that by performing an overall best-fit (instead of using 
the "data" for on-shell momenta k < 0.1 fm _1 ) a fairly good description of the phase- 
shifts can be obtained up to Slab ~ 200 MeV, but that would prevent us from making a 
proper error scaling analysis for predicted phase-shifts. Furthermore, as a consequence of 
the Wigner causality bound [79] we find that by taking a subtraction scale /i larger than 
~ 1.1 fm" 1 it becomes impossible to fit the inverse on-shell K- matrix evaluated from 
the Nijmegen-II potential to within an acceptable accuracy while keeping the NNLO 
driving term hermitian (i.e., with the renormalized strengths Cj(/i) taking only 
real values), such that the unitarity of the scattering amplitude is preserved. This is 
illustrated by the results obtained for /i = 1.2 fm -1 , which cannot reproduce well the 
Nijmegen-II phase-shifts even at very low energies, and so clearly indicate the breakdown 
of the fitting procedure. Several works have discussed the implications of the Wigner 
causality bound [2TI EH [531 EOl EU [821 [831 EH [85] and the so called low-energy theorems 
[221 l86l [87] to the renormalization problem, both in the context of pionless EFT and 
ChEFT. The results described in the works cited above give strong support to the 
prescription advocated by Lepage [8] that in order to get a consistent cutoff EFT for 
the NN interactions the cutoff should not be taken much larger than the relevant hard 
scale in the theory, e.g. the pion mass m n ~ 140 MeV in the case of pionless EFT 
and the chiral symmetry breaking scale A x ~ 1 GeV in the case of ChEFT. The most 
efficient choice is to take the cutoff of the order of the relevant hard scale [2j El [2H [22] . 

EFT's are essentially constructed as systematic low-energy expansions in powers 
of the ratio (Q/Ahard), where Q and Ahard stand respectively for the relevant soft (low- 
momentum) and hard (high-momentum) scales in the theory, which provide a valid 
description of phenomena at momentum scales below Ahard- In the renormalization 
of a cutoff EFT for the NN interactions, an UV regularizing momentum cutoff scale 
A is introduced in order to remove high-momentum degrees of freedom which can 
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probe the unknown short- distance dynamics, thus playing the role of a resolution 
scale. The contributions from low- momentum/long-distance degrees of freedom are 
included explicitly in the calculation of observables through interactions known from 
the underlying theory (e.g. pion-exchange interactions). The contributions from the 
excluded high-momentum/short-distance degrees of freedom are included implicitly 
through parametrized cutoff-dependent contact interactions (counterterms) whose 
coefficients are fixed by fitting low-energy scattering data, systematically removing cutoff 
dependence in the observables. Lepage's prescription is based on the view that by 
taking the momentum cutoff scale A much larger than A^ard (i-e., beyond the range of 
validity of the EFT) one is certainly incorporating contributions from incorrect high- 
momentum/short-distance dynamics which can lead to pathologies, such as the violation 
of the Wigner causality bound, and even the breakdown of the EFT systematics j22j EH] . 

The Wigner causality bound is a general result, originally derived by Wigner [79] 
assuming only the physical principles of causality and unitarity, which shows that for a 
hermitian potential that vanishes beyond some range R there is a lower bound on the 
derivative of the phase-shifts S(k) with respect to the on-shell momentum k, 



An alternative derivation was presented later by Phillips and Cohen |80j, who have 
shown that Wigner 's bound yields a constraint on the effective range r e , given by 



where a is the scattering length. Moreover, it was shown that this constraint still applies 
when the potential is not identically zero beyond the range R, but fall off sufficiently 
fast for the wave-function to approach its asymptotic solution rapidly enough. 

Eq. (123 p shows that in the limit of a zero-range interaction (R — > 0) the effective 
range r e cannot be positive. It also shows that there is a minimum range -R m i n for which 
the interaction can reproduce a given scattering length a and effective range r e . In 
the case of NN scattering in the 1 Sq channel, for which the experimental values of the 
scattering length and the effective range are respectively given by a = —23.7 fm and 
r e = 2.77 fm, this minimum range is -R m i n = 1-3 fm. Furthermore, Scaldeferri et. al. 
[8Tj have shown that this constraint is a general feature of contact interactions which 
still holds even when long-range pion-exchange interactions are explicitly included. In 
particular, it was shown that an absolute lower bound -R m i n = 1.1 fm is obtained for 
the AW scattering in the 1 5'o channel when both the OPE interaction and an arbitrary 
short-range contact interaction are included. 

By identifying the range R with the momentum cutoff scale A introduced in the 
renormalization of cutoff EFT for the AW interactions (i.e., R ~ 1/A), Wigner's bound 
implies both for pionless EFT and ChEFT that in order to obtain a description of phase- 
shifts which agrees with experimental low-energy AW scattering data, the momentum 
cutoff scale cannot be removed by taking the limit A — > oo while maintaining causality 



d5(k) 
dk 



>-R+-y sm(25(k) + 2kR) . 



(22) 




(23) 
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and unitarity, as required for a consistent EFT framework J27J E2J E2J El]- Indeed, 
by investigating the renormalizability of chiral NN interactions in the 1 5'o channel at 
NLO and NNLO in the WPC scheme within the framework of cutoff renormalization, 
Entem et.al. [27] have found that there is a maximum value A max for the momentum 
cutoff scale A above which one cannot fix the strengths of the contact interactions 
C (A) and C 2 (A) by fitting the experimental values of both the scattering length a 
and the effective range r e while keeping the renormalized potential hermitian. For 
A > A max ~ 1.8 fm _1 in the case of pionless EFT and A > A max ~ 2.5 fm _1 in the case 
of ChEFT, the strengths Co (A) and C2(A) diverge before taking complex values. Thus, 
in both cases the renormalized potential must become non- Hermitian in order to match 
the renormalization conditions, such that the unitarity of the scattering amplitude is 
violated (nevertheless, the corresponding phase-shifts remain real). These results are 
similar to those we have obtained in our calculations for chiral AW interactions in the 1 Sq 
channel at NLO and NNLO within the framework of the SKM approach. The rather 
smaller value for the maximum subtraction scale /i max ~ 1.1 fm _1 we have found in the 
SKM approach as compared to A max ~ 2.5 fm -1 found in the cutoff renormalization 
scheme is a consequence of the highly non-trivial energy, momentum and subtraction 
scale dependence of the scattering amplitude K^> obtained from the solution of the 
n-fold subtracted LS equation, Eq. (fl3|) . which yields a very different running of the 
renormalized strengths of the contact interactions C(yu) with the subtraction scale \x as 
determined by the renormalization conditions. 

It is important to emphasize that both in the cutoff renormalization scheme 
and the SKM approach, Wigner's bound arises from the non-linear structure of the 
renormalization conditions relating the strengths of the contact interactions C to the 
ERE parameters (a, r e , v 2 , etc) used as physical input [821 [831 El]- Therefore, the 
constraint on the maximum value of the regularizing scale for which one can fit the 
experimental values of the ERE parameters (beyond zeroth order) with a hermitian 
potential (i.e., with the strengths of the contact interactions Cj taking only real values) 
applies regardless of the particular procedure used to implement the fitting. 

4. Renormalization group invariance in the SKM approach. 

As pointed before, the multiple subtractions performed in the SKM procedure introduce 
a form factor in the kernel of the LS equation which acts like a regularizing function, 
such that the subtraction scale /i ends up playing a role similar to that of a smooth 
momentum cutoff scale. But the subtraction scale /i is arbitrary, and so the scattering 
observables calculated from the solution of the subtracted LS equation for the scattering 
amplitude should not depend on its particular choice. By requiring the fully off-shell 
if-matrix with n subtractions to be invariant under the change of the subtraction scale 
/i, a renormalization group (RG) equation can be derived for the driving term V^ n \E) 
in the form of a non-relativistic Callan-Symanzik (NRCS) flow equation [57], which is 




Figure 3. (Color on-line) Evolution through the NRCS equation of the 1 5'o channel 
driving term at NNLO in the MPC scheme for several values of Slab- Top 
panels: diagonal matrix elements; Bottom panels: off-diagonal matrix-elements. 



given in operator form by 

^-^'(^^(^ P4) 

with the boundary condition V^\^p_ = imposed at some reference subtraction 
scale ft where the renormalized strengths of the contact interactions Cj(/i) are fixed to 
fit low-energy observables used as physical input. 

The NRCS flow equation for the driving term is similar to the RG equation for 
the NN potential derived by Birse in the context of cutoff EFT [IHJ H9j [36] , which is also 
obtained from the invariance of the off-shell .fT-matrix. Another similar RG equation is 
that derived by Bogner et al. in the V[ ow _k approach (HHl IHSl ED] , which is based on the 
invariance of the half-on-shell T-matrix such that it also involves the iteration of the 
scattering amplitude. 

In order to explicitly demonstrate the renormalization group invariance in the SKM 
approach, we consider the evolution through the NRCS flow equation of the 1 5o channel 
driving term at NNLO in the MPC scheme. In a partial-wave relative momentum 
space basis, the NRCS flow equation for the matrix-elements of the driving term 
in uncoupled channels is given by 



dVi n \p,p';k 



2\ 



(fx 2 + k 

n 



2\n-l 



f°° 2 

<V n'o dqq f(/i 2 + g 2 ) n+1 . 

xV^^q-k^V^^p'-k 2 ) . (25) 

We solve Eq. f[2"5"j) numerically for n = 4, obtaining an exact (non-perturbative) 
solution for the evolved 1 Sq channel driving term V^> . The relative momentum space is 
discretized on a grid of 200 gaussian integration points, leading to a system of 200 x 200 
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Figure 4. (Color on-line) Phase-shifts in the l So channel obtained with the driving 
term at NNLO in the MPC scheme evolved through the NRCS equation. 

For comparison, we also show the phase-shifts obtained with the driving term 
computed at the reference scale p, (solid line). 



non-linear first-order coupled differential equations which is solved using an adaptative 
fifth-order Runge-Kutta algorithm. In Fig. [3] we show the evolution of the diagonal 
matrix-elements (top panels) and the off-diagonal matrix-elements (bottom panels) of 
the 1 5'o channel driving term from a reference subtraction scale fx = 1.03 fm -1 
to /i = 0.83, 0.93 and 1.2 fm -1 for several values of -Elab- One should note that the 
solution of the NRCS flow equation leads to a non-trivial evolution of the driving term 

with the sliding subtraction scale /i. The change in the form of the driving term 
is not only due to the running of the renormalized strengths of the contact interactions 
Ci(fi) but also due to the new operators that are generated by the NRCS flow. 

As shown in Fig. HI the evolution of the 1 Sq channel driving term Vff^ through 
the NRCS equation ensures that the phase-shifts calculated from the solution of the 
LS equation for the 4-fold subtracted .fT-matrix remain invariant (except for relative 
differences smaller than 10~ 12 due to numerical errors). Remarkably, this result still 
holds even when the driving term is evolved to subtraction scales fx which are beyond 
the limit imposed by the Wigner causality bound (such as fx = 1.2 fm -1 ). This can 
be understood in the following way. As pointed before, the input physical information 
(e.g. the experimental values of the ERE parameters a, r e , t> 2 , etc) is encoded in the 
initial driving term V^ n \ which is determined by fixing the renormalized strengths Ci(fi) 
of the contact interactions at the reference subtraction scale fx to fit data and sets the 
boundary condition for the NRCS evolution. Thus, once the constraint imposed by 
Wigner's bound is fulfilled at the boundary condition, provided one chooses a reference 
scale fx < Umax (which in this case is ~ 1.1 fm" 1 ), the NRCS flow equation, Eq. ( 125]) . 
can be solved to evolve the driving term to any subtraction scale fx while keeping it 
hermitian. These results clearly show that the sliding subtraction scale fx vanishes as 
a physical parameter, the only relevant scale being the reference scale fx where the 
boundary condition for the NRCS evolution is determined. 
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5. Summary and conclusions. 

The subtracted kernel method (SKM) approach presented in this work provides 
a powerful renormalization group invariant method to renormalize singular two- 
body interactions. The iterative procedure involving recursive multiple subtractions 
performed in the kernel of the scattering equation allows for a systematic treatment 
of chiral effective field theory (ChEFT) nucleon-nucleon (NN) potentials up to higher- 
orders in the chiral expansion, which include pion-exchange and contact interactions 
that present increasingly strong singularities at short-distances. This also makes the 
SKM approach a convenient tool to implement power-counting schemes in which higher- 
order contact interactions are promoted to lower-order, since ultraviolet (UV) power 
divergences of any order can be properly handled by performing a number of subtractions 
enough to render a finite amplitude. A disadvantage of the SKM formalism is that the 
computational load increases as more subtractions are performed: in order to compute 
the n-fold subtracted amplitude, the iterative procedure requires n matrix inversions. 

We have shown that a modified power counting scheme (MPC) based on the 
promotion of contact interactions, such that at each order in the chiral expansion one 
new contact interaction is included in the effective NN potential, yields a systematic 
order- by-order power-law improvement in the 1 So channel phase-shifts calculated up to 
next-to- next-to-leading-order (NNLO), that is lacking in Weinberg's power counting 
scheme (WPC). Our results show that the scaling of the relative errors at low- 
energies is nearly dominated by the tuning of the renormalized strengths of the contact 
interactions, and so suggest that the short-range contact interactions are relatively more 
important than the long-range pion-exchange interactions. An essential ingredient in our 
calculations is the procedure employed to fix the renormalized strengths of the contact 
interactions, based on the fitting of "data" generated from the Nijmegen II potential 
for a spread of very low on-shell momenta, which allows for a clear analysis of the error 
scaling of predicted phase-shifts, and hence of the power counting, using "Lepage plots" . 
Such a procedure is equivalent to a fitting of the experimental values of the effective 
range expansion (ERE) parameters (a, r e , v 2 , etc) to a given order in k 2 , where k is the 
on-shell momentum in the center-of-mass frame. 

One should note that we are not analyzing the power counting at the level of the 
chiral expansion for the NN potential, but "a posteriori" at the level of the observables, 
i.e. after iterating the potential in the subtracted Lippmann-Schwinger (LS) equation. 
The lack of power- law improvement in the 1 Sq channel phase-shifts we have observed at 
NNLO in the WPC scheme shows that the corresponding subtracted amplitude do not 
follow the power counting of the chiral NN potential, which is the motivation for our 
MPC scheme. As pointed out by several works [H [22j [371 [381 E3 E0 EE], a consequence 
of the singular nature of the chiral NN interactions at short-distances is that the full 
iteration of the potential in the scattering equation can change the scaling behavior of 
each of the interactions included, and hence their relative importance to the amplitude, 
thus leading to modifications of the power counting at the level of the observables. 



Power counting and renormalization group invariance in the SKM 



16 



We have also demonstrated by explicit numerical calculations for the scattering 
of two-nucleons in the 1 Sq channel that the SKM procedure applied to the chiral NN 
interactions up to NNLO is renormalization group invariant under the change of the 
subtraction scale. Once the renormalized strengths of the contact interactions are fixed 
at a reference scale to fit low-energy scattering observables, the subtraction scale can 
be changed by evolving the driving term of the subtracted LS equation through a 
non-relativistic Callan-Symanzik (NRCS) flow equation, such that the results for the 
calculated phase-shifts remain invariant. Moreover, the sliding subtraction scale to 
which the driving term can be evolved is not constrained by the Wigner causality bound. 
In this way, the sliding subtraction scale vanishes as a physical parameter. The relevant 
scale parameter left in the theory is the reference scale where the boundary condition of 
the NRCS flow equation is determined through the input of physical information, which 
can then be regarded as a renormalization scale. 

It is important to emphasize that in order to properly assess the effectiveness of 
the MPC scheme we have described in this work, as well as the renormalization group 
invariance in the SKM approach, it is necessary to carry out the calculations up to 
next-to-next-to-next-to-leading-order (N 3 LO) and perform a comprehensive analysis 
for all partial-wave channels relevant at low-energies (e.g. for total angular momentum 
j < 5). In particular, the calculations for the lower partial- waves would be essential to 
address the consequences of our results for applications of ChEFT to few- and many- 
nucleon systems (e.g. light-nuclei, and nuclear and neutron matter). If we do find, 
for all relevant partial- wave channels and up to N 3 LO, that a systematic order-by- 
order power-law improvement can be achieved through the MPC scheme and the NN 
phase-shifts remain invariant under the change of the subtraction scale through the 
NRCS flow equation, then it would be reasonable to expect that similar results may 
also be obtained when applying the SKM approach to few- and many-nucleon systems. 
Of course this must be verified through explicit calculations, which would require the 
extension of the SKM approach to solve dynamical equations pertinent to nuclear few- 
and many-body problems (e.g. the Faddeev-Yakubovsky equations [91] and the Bethe- 
Brueckner-Goldstone equations [92]). Nevertheless, due to its sensitivity to the short- 
range contributions to the chiral NN potential [27] [51] , the 1 5'o channel provides a good 
starting point to investigate the interplay between the pion-exchange and the contact 
interactions in the non-perturbative renormalization of NN interactions in ChEFT. 

In forthcoming works, we will extend the calculations presented here to other 
partial- wave channels and up to N 3 LO, aiming to perform a systematic analysis of 
the power counting and the renormalization group invariance in the SKM approach 
for chiral NN interactions. In particular, we will investigate in detail the impact 
of the long-range pion-exchange interactions and their interplay with the short-range 
contact interactions. We also intend to compare the evolution of chiral NN interactions 
through the NRCS equation in the SKM approach with the evolution through Birse's 
renormalization group (RG) equation in cutoff EFT, both in a perturbative and in a 
non-perturbative renormalization framework. 
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Furthermore, we want to investigate the possibility of generalizing the SKM 
approach to renormalize ChEFT potentials including three- and four-nucleon 
interactions. We believe that this generalization can be worked out by using techniques 
similar to those developed in the context of EFT for three-body systems with contact 
interactions [931 EH ESI EE], which allow to obtain renormalized amplitudes by 
performing subtractions in the three-body scattering equations. By accomplishing such 
a generalization, we expect to be able to apply the SKM approach in sensible calculations 
of few- and many-nucleon systems. 
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